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Summary 

This report documents the results of a study conducted to compute the inviscid longitudinal aerodynamic* 
characteristics of a simplified X-33 configuration. This study was conducted in support of the industry led 
X-33 program. The major components of the X-33 vehicle, namely the body, the canted fin, the vertical fin, 
and the bodyflap, were simulated in the CFD model. The rear-ward facing surfaces at the base including the 
aerospike engine surfaces were not simulated. The FELISA software package was used for this study. This 
software consists of an unstructured surface and volume grid generator and two inviscid flow solvers one for 
transonic flows and the other for hypersonic flows. The hypersonic flow solver with perfect gas, equilibrium 
air, and CF 4 gas options was used in the present study. Computations were made for Mach 4.96, 6.0. 
and 10.0 with perfect gas air option, and for Mach 10 with equilibrium air option with flow condition of a 
typical point on the X-33 flight trajectory. Computations were* also made with CF 4 gas option at Mach 6.0 
to simulate the CF 4 tunnel flow condition. An angle of attack range of 12 to 48 deg. was covered. The 
CFD results were compared with available 1 wind tunnel data. Comparison was good at low angles of attack; 
at higher angles of attack (beyond 25 deg.) some differences were found in the pitching moment. These 
differences progressively increased with increase in angle of attack, and are attributed to the viscous effects. 
However, the computed results showed the trends exhibited by the wind tunnel data. 


Nomenclature 


C A 

C\ 

c m 

c„ 


Fa /( Qoc $ref), Axial force coefficient 
Fn/( q 0 c S re f), Normal force coefficient 
M y/(<7oc $ref kef), Pitching moment coefficient 
(p - p 0 o)/<7oo 5 Pressure coefficient 


Fa 

Fn 

l re f 

My 

V 

Pa c 
Qoo 

Sref 
x, y, z 

(X 


Axial force 
Normal force 

Reference length, ( = 758.4sq.in.) 
Pitching moment 
Freestream Mach number 
Static pressure 
Freestream static pressure 
Freestream dynamic pressure 
Reference area, ( = 231,520.0sq.in.) 
Cartesian co-ordinates of a given point 
Angle of attack, deg. 


Geometry 

The present computational study was done on the Lockheed X-33 F-loft, Rev. F configuration. The primary 
external components of this vehicle are the body, the canted fin, the vertical fin, the bodyflap, and the 
aerospike engines. All of these components, except the aerospike engine surfaces, were simulated in the CFD 
model. In order to preclude separated flow regions, mainly over the rearward facing areas like the base of 
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Figure 1: The X-33 F-loft Rev. F vehicle. 


the vehicle, the flow over the base was not simulated. This was accomplished by extending the surfaces of 
the canted fin, vertical fin, and the of body downstream of the trailing edge so that they ended on a single 
surface. These extension surfaces provide a cover for the base. In the CFD computations, the flow was 
forced to be tangential to these extension surfaces. However, while computing the aerodynamic loads on the 
vehicle, these surfaces were ignored. 

A sketch of the X-33 model used in the present CFD study is shown in Fig. 1. The axis system used 
for this study is also shown in the figure. The origin is at the nose of the vehicle with the z-axis along the 
body axis pointing upstream, y-axis perpendicular to the symmetry plane, and the x-axis in the symmetry 
plane perpendicular to- the z-axis. The length of the body (from nose to trailing edge of the bodyflap) is 
approximately 831.1 in. and the tip-to-tip distance of the canted fins is 920.75 in. The canted fins make 
an angle of 72 deg. to the plane of symmetry. There are two vertical fins separated from each other by 
200 in. For the present computations, the control surfaces, namely the rudder and the bodyflaps were 
in the undeflected position. Since the vehicle has a plane of symmetry, only one half of the vehicle was 
simulated in the computational model. The reference quantities used for reducing the aerodynamic loads to 
the non-dimensional form are as follows: 


Reference area 

Reference length 

Pitching moment reference point 


231,520.0 sq.in. 

758.4 in. 

(0.0, 0.0, -500.544) in. 


Grid ID 

No. of Tets 

No. of Points 

No. of Triangles 

No. of Surface Points 

X33A 

7,955,084 

1.337,224 

242,850 

121,427 

X33C 

6,506,154 

1,094,694 

203,160 

101,582 

X33D 

7,069,961 

1.189,736 

220.066 

110,035 


Table 1: Properties of grids used in the present computations. 


The Felisa Software 

All the computations of the present study were done using the FELISA unstructured grid software. This 
software package consists of a set of computer codes for the simulation of three-dimensional steady inviscid 
flows using unstructured tetrahedral grids. Surface triangulation and discretization of the computational 
domain using tetrahedral elements is done by two separate codes. Two flow solvers are available -one 
applic able for transonic flows, and the other for hypersonic flows with an option for perfect gas air, equilibrium 
air, or CF 4 gases. The hypersonic flow solver was used for all the computations in this study. Post-processors 
like the aerodynamic analysis routine are part of the software package. More information on FELISA may 
be found in [1], A description of the hypersonic flow 7 solver may be* found in [2]. 

Grids Used in the Present Study 

Stai ting with a structured surface grid in the form of a PLOT3D file, and using the software GridTool [3], 
a set. of FELISA data files was generated. This includes the body surface and the computational domain 
definition file, and a file that specifies the grid spacing. The grid spacings file was modified manually in 
ordei to obtain the desired grid spacing near the nose of the vehicle. Using these data, files unstructured 
tetrahedral grids were generated with the FELISA volume grid generator. Three separate unstructured 
surface and volume grids were generated to suit the freestream conditions. The properties of these grids are 
shown in Table 1 above. All these grids were generated on an SGI ONYX located in the Aerothermodvnamics 
Branch (AB), NASA Langley. 

A minimum grid spacing was 0.4 inch was used near the nose of the body. A plot of a portion of the grid 
on the symmetry plane near the nose is shown in Fig. 2. The model surface triangulat.ion for grid X33A is 
shown in Fig. 3. 


Flow Conditions 

The freestream conditions and the range of angles of attack covered in the present study are shown in 
Table 2. All perfect gas air computation (cases 1 12) were done with 7 = 1.4. Freestream conditions for 
the equilibrium air computations (cases 13-15) corresponded to the conditions for a typical point on the 
X-33 flight trajectory. Some wind tunnel tests had been conducted on an X-33 model in the Langley 31-inch 
Mach 10 tunnel and also in the Langley Mach 6 CF 4 tunnel. In order to simulate the CF 4 tunnel tests, a set 
of computations was done (cases 18-21) with the CF 4 gas options and freestream conditions corresponding 
to the conditions in the CF 4 tunnel. The freestream conditions for computations at. Mach 9.02 correspond 
to atmospheric flight at a. nominal Mach number of 9 at an altitude of 133,000 ft.., and the conditions for 
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Figure 3: Surface triangulation for X33A grid. 
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Case No. 

Mach No. 

Gas Model 

Density 

(Kg/m 3 ) 

Temperature 

(Kelvin) 

Velocity 

(m/sec) 

1-4 

10.0 

Perfect Gas 

NA 

NA 

NA 

5 8 

6.00 

Perfect Gas 

NA 

NA 

NA 

9-12 

4.96 

Perfect Gas 

NA 

NA 

NA 

13-15 

10.0 

Equil. Air 

4.3120E-04 

259.5 

3224 

16 

9.02 

Equil. Air 

3.6867E-03 

251.9 

2863 

17 

12.02 

Equil. Air 

1.0593E-03 

270.7 

3958 

18-21 

6.02 

CF 4 gas 

1.7700E-02 

170.0 

850 


Table 2: Freest ream conditions . 


computations at Mach 12.02 correspond to atmospheric flight at a nominal Macli number of 12 at an altitude 
of 163,200 ft. These two cases correspond to two critical points in the flight trajectory of the X-33 vehicle. 

Flow Solution 

The flow 7 solutions were computed on SGI Origin 2000 class parallel computers having 64 processors sitting 
on top of 16G of shared memory. The grids were partitioned so that the problem would run on eight 
processors. The FELISA hypersonic flow solver with the appropriate gas option w 7 as used for these flow 7 
computations. Each solution was started with the low-order option, and after a few 7 hundred iterations, 
the higher-order option was turned on, and the solution was run to convergence. After every 20 iteration, 
the surface pressures were integrated, and the aerodynamic loads, namely the normal and the axial forces, 
and the pitching moment acting on the model w 7 ere computed. The flow solution was assumed converged 
when these integrated loads reached steady values. Each solution required 32-40 hours of CPU time. The 
computed flow' solutions were post-processed on the AB ONYX computer to obtain the aerodynamic loads. 
These loads w 7 ere non-dimensionalized in the conventional manner, and the aerodynamic coefficients (Cyv, 
Caj and C m ) were obtained. 

It should be recalled at this point that all the present computations are inviscid; hence the boundary 
layer and skin friction are absent. Absence of the skin friction leads to lower axial forces. Further, since the 
boundary layer is absent, the effects of boundary layer separation on the aerodynamic coefficients are not 
simulated. On the X-33, the flow 7 would be fully separated over the base, and possibly over the upper surface 
of the bodyflap. These surfaces are not simulated in the present CFD model. These factors could introduce 
inaccuracies in the computed axial and normal forces, and as well as in the pitching moment. 


Results and Discussion 

The results of the present study are summarized in Table 3, and also shown graphically in Figures 4, 5, and 
6. Tunnel test data in the form of curve fits for Cat, Ca, and C rn were available for the cases listed in the 
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table 4. For the Mach 6 air case, the data at 4 million Reynolds number was chosen for comparison with the 
computed results. 

The computed CVv and C m for Mach numbers 4.96, G, and 10 in perfect gas air are shown plotted Vs. a 
in Figure 4. Smooth curves have been drawn through the data points. The C N Vs. a curves for the three 
cases shown in Fig. 4(a) have nearly the same slope. For a given a, the C\ values decrease as the Mach 
number is increased from 4.96 to 10.0. The pitching moment data for the three Mach numbers is shown in 
Fig. 4(b). For Mach 4.96 the C m Vs. o has a large negative slope for a up to 35 deg.; beyond this angle 
of attack the slope becomes less negative. This indicates that the model is stable up to o=35 deg., beyond 
which it becomes less stable. For Mach 10. the trend is different . The slope of C m Vs. a curve is nearly 
zero up to o=25 deg.; beyond this angle of attack the slope is negative. This indic ates that the model is 
neutrally stable up to a =25 deg., beyond this angle of attack the vehicle becomes progressively more stable. 
The C m Vs. a for Mach 6 lies in between the curves for Mach 4.96 and 10.0. 

Computed Cj\ and C m values for Mach 10 perfect gas air arid equilibrium air are shown in Fig. 5. 
Tunnel data for Mach 10 are also plotted in this figure for comparison. Figure 5(a) shows that there is little 
difference between the computed Ck values for perfect gas air and the equilibrium air assumption. It should 
be recalled here that the freestream conditions for these equilibrium flow computations corresponded to those 
for a typical point on tin* X-33 flight trajectory. The total temperature 1 under these conditions is 3040K. At 
this high temperature there would be considerable changes in the chemical composition and properties of air 
(see e.g. [4]). These changes do not seem to affect the normal force coefficient. The total temperature for 
the freestream conditions in the Mach 10 tunnel test section is only 980K. At this temperature there is very 
little change in the chemical composition of air; hence these tunnel data should be compared with perfect 
gas computations. 

Figure 5(b) shows that there are some differences between computed C m values for perfect gas and 
equilibrium air cases. These differences are due to the differences between the properties of perfect gas air 
and equilibrium air (at a total temperature of 3040K). However, there are significant differences between the 
computed and the wind tunnel C m data. The differences are, although small at low angles of attack (0.001 
up to o— 30 deg.), increase as the angle of attack is increased (0.004 at o=48 deg.). This trend is possibly 
due to the viscous effects. On the wind tunnel model the flow on the upper surface of the canted fin would 
separate at high angles of attack due to the presence of the trailing edge shock. This would increase the 
pressures over the aft part of the canted fin, and, due to the large 1 moment arm, result in positive pitching 
moment. 

The computed C N and C m for Mach 6 perfect gas air and CF 4 gas are shown in Figs. 6(a) and (b) 
respectively, along with the test data from Langley 20-Inch Mach 6 Air tunnel and Langley Mach 6 CF 4 
tunnel. The freestream conditions for the CF 4 gas computations corresponded to those in the Langley Mach 
6 CF 4 tunnel test section with a total temperature of 980K. Figure 6(a) shows that the assumption of CF 4 
gas at Mach 6 does not affect the computed normal force. Further, there is good agreement, between the 
computations and the CF 4 tunnel data. However, there are differences between the computed results and 
the air tunnel data. In the past Cat values computed using FELISA had compared well with tunnel data; 
see for example [5]. The differences noticed in the present case could not be explained. 

The computed C m values for perfect gas air and CF 4 gas cases for Mach 6 are shown in Fig. 6(b). It may 
be observed that there are differences between the C m values for the two cases. This is primarily due to the 
fact that air and CF 4 gas have different thermodynamic properties. For example, the density ratios across 
a normal shock at Mach 6 in CF 4 gas is about 10 compared to a value of about 5.3 in air. The pressure 
distribution over the vehicle for the CF 4 gas and in perfect gas air would be different. The tunnel test data 
is also shown in Fig. 6(b). It may be noticed from this figure that the computed C m for CF 4 gas compares 
with the tunnel data up to about a=30 deg. Beyond this angle, the tunnel data exhibits less nose down C m 
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Case No. 

Mach No. 

a, deg. 

Gas Model 

C N 

Cm 

Cm 

1 

10.0 

24 

P.G. Air 

0.4534 

0.1033 

0.0105 

2 

10.0 

36 

P.G. Air 

0.8922 

0.1014 

0.0037 

3 

10.0 

48 

P.G. Air 

1.3245 

0.0919 

-0.0033 

4 

10.0 

12 

P.G. Air 

0.1232 

0.1048 

0.0111 

5 

6.00 

24 

P.G. Air 

0.4898 

0.1030 

0.0089 

6 

6.00 

36 

P.G. Air 

0.9252 

0.0998 

0.0032 

7 

G.OO 

48 

P.G. Air 

1.3401 

0.0909 

-0.0002 

8 

6.00 

12 

P.G. Air 

0.1434 

0.1065 

0.0118 

9 

4.96 

24 

P.G. Air 

0.5146 

0.1025 

0.0082 

10 

4.96 

36 

P.G. Air 

0.9481 

0.0979 

0.0030 

11 

4.96 

48 

P.G. Air 

1.3555 

0.0881 

0.0012 

12 

4.96 

12 

P.G. Air 

0.1579 

0.1071 

0.0123 

13 

10.0 

36 

Equil. Air 

0.8853 

0.1079 

0.0048 

14 

10.0 

48 

Equil. Air 

1.3286 

0.1040 

-0.0041 

15 

10.0 

24 

Equil. Air 

0.4484 

0.1070 

0.0113 

16 

9.02 

20 

Equil. Air 

0.3279 

0.1062 

0.0116 

17 

12.02 

35 

Equil. Air 

0.8318 

0.1092 

0.0060 

18 

6.02 

12 

CF 4 Gas 

0.1458 

0.1090 

0.0118 

19 

6.02 

24 

CF 4 Gas 

0.4763 

0.1065 

0.0116 

20 

6.02 

36 

CF a Gas 

0.9044 

0.1057 

0.0059 

21 

6.02 

48 

CF 4 Gas 

1.3526 

0.0991 

-0.0039 


Table 3: Summary of computed aerodynamic coefficients for X-33. 


Tunnel 

Mach No. 

Reynolds No. 

a Range 

Mount 

Langley 20-Inch Mach 6 Air 

6.0 

0.5E+06 

-4 to 48 deg. 

Sting 

Langley 20- Inch Mach 6 Air 

6.0 

2.0E+06 

-4 to 48 deg. 

Sting 

Langley 20-Inch Mach 6 Air 

6.0 

4.0E+06 

-4 to 48 deg. 

Sting 

Langley CF A 

6.0 

.0.4E+06 

-4 to 48 deg. 

Sting 

Langley 31-Inch Air 

10.0 

2.0E+06 

5 to 50 deg. 

Blade 


Table 4: Available wind tunnel data for X-33. 
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(a) C N Vs. o 



a, deg. 


(b) Cm Vs. a 


Figure 4: Computed CV and C rn for X-33 in perfect gas air at Mach 4.96, 6.0, and 10.0. 
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a 

Body 

Canted Fin 

Bodyflap 

(deg.) 

C N 

C A 

c m 

c N 

C A 

Cm 

c N 

C A 

c 

^ m 

12 

0.1221 

0.0834 

0.0137 

0.1221 

0.0834 

0.0137 

0.1221 

0.0834 

0.0137 

24 

0.3748 

0.0931 

0.0325 

0.3748 

0.0931 

0.0325 

0.3748 

0.0931 

0.0325 

36 

0.7182 

0.1023 

0.0518 

0.7182 

0.1023 

0.0518 

0.1067 

-0.0025 

-0.0242 

48 

1.0691 

0.1050 

0.0678 

1.0691 

0.1050 

0.0678 

0.1619 

-0.0147 

-0.0380 


Table 5: Computed aerodynamic data for body, canted fin, and bodyflap, Mach 10, perfect gas air. 


a 

Body 

Canted Fin 

Bodyflap 

(deg.) 

C N 

C A 

c m 

C'N 

Ca 

c,„ 

Cn 

Ca 

c m 

24 

0.3715 

0.0965 

0.0329 

0.0427 

0.0094 

-0.0095 

0.0342 

0.0010 

-0.0122 

36 

0.7098 

0.1085 

0.0535 

0.1058 

-0.0018 

-0.0240 

0.0695 

0.0012 

-0.0247 

48 

1.0608 

0.1170 

0.0706 

0.1655 

-0.0138 

-0.0386 

0.1021 

0.0008 

-0.0361 


Table 6: Computed aerodynamic data for body, canted fin, and body flap, Mach 10. equilibrium air 


. For the perfect gas air case, the comparison is not as good. The difference increases progressively with 
increase in a for zero at a=12 deg. to 0.006 at, o=48 deg. 

In the Mach 6 and Mach 10 cases, the tunnel C m values are invariably higher than the computed values 
at angles of attack 30 deg. and beyond. As noted earlier, this is possibly due to the effects of flow separation 
on the canted fin. The present computations are inviscid, as such there is no flow separation any surface, 
including the canted fins. 

The contributions to Cs and C m due to the body, canted fin, bodyflap, and the vertical fin were computed 
separately for all the cases. In all the cases the contributions from the vertical fin to C s and C in were found 
to be negligible. The contribution from the body, canted fin, and the bodyflap for Mach 10 perfect gas and 
Mach 10 equilibrium air cases are listed in Tables 5 and 6, and are also shown plotted in Figures 7 and 
8. It may be noticed from these figures that body contributes the most to both normal force and pitching 
moment. The contributions of the canted fin and the bodyflap are of comparable magnitude. 

Similarly, the results for Mach 6 perfect gas air and Mach 6 CF 4 gas cases are listed in Tables 7 and 8, 
and art' also shown plotted in Figure 9 and 10. 


Conclusion 

Inviscid longitudinal aerodynamic characteristics of the Lockheed X-33 f-loft, Rev. F vehicle were computed 
for Mach 4.96, 6.0, and 10.0 over an angle of attack range of 4 -48 deg. using the FELISA unstructured grid 
software package. The bodyflap and the rudder were in undeflected position. The canted fin, the vertical fin 
and the body surfaces were extended downstream in the computational model in order to preclude regions 
of separated flows. Computations were made for Mach 10 with the perfect gas air and equilibrium air 
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(a) C\ Vs. a 



(b) Cm Vs. a 


Figure G: Comparison of computed Cn and C m for X-33 in perfect gas air and CF A gas at Mach 6 
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(a) Op/ Vs. a 



(b) Cm Vs. Q 


Figure 7: Computed C/v and C m due to the body, canted fin, and bodyflap, and the complete vehicle for 
M— 10, perfect gavS air. 
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(b) Cm Vs. a 










o 

Body 

Canted Fin 

Bodyflap 

(deg.) 

Cj\ 

C A 

c m 

C N 

c A 

c m 

C N 

C A 

Cn, 

12 

0.1439 

0.0848 

0.0140 

-0.0115 

0.0202 

0.0017 

0.0113 

0.0008 

-0.0041 

24 

0.4018 

0.0940 

0.0336 

0.0510 

0.0078 

-0.01 1G 

0.036G 

0.0010 

-0.0130 

36 

0.7405 

0.1029 

0.0539 

0.1159 

-0.0042 

-0.0264 

0.0682 

0.0011 

-0.0242 

48 

1.0791 

0.1065 

0.0720 

0.1670 

-0.0162 

-0.0391 

0.0936 

0.0001 

-0.0330 


Table 7: Computed aerodynamic data for body, canted fin, and bodyflap, Mach G. perfect gas air. 


a 

Body 

Canted Fin 

Bodyflap 

(deg.) 

C N 

C A 

c,„ 

Cn 

C A 

c 

in 

C'n 

C A 

c 

171 

12 

0.1 45G 

0.0873 

0.0143 

-0.0105 

0.0202 

0.0013 

0.0110 

0.0008 

-0.0039 

24 

0.3977 

0.09G0 

0.0339 

0.0441 

0.0093 

-0.0100 

0.0341 

0.0010 

-0.0122 

3G 

0.7280 

0.1058 

0.0550 

0.1062 

-0.0015 - 

-0.0242 

0.0697 

0.0013 

-0.0247 

48 

1.0796 

0.1118 

0.0720 

0.1695 

-0.0138 - 

-0.0393 

0.1028 

0.0008 

-0.0309 


Table 8: Computed aerodynamic data for body, canted fin, and bodyflap, Mach 6. CF 4 gas. 


assumptions, and for Mach G with the perfect gas air and CF 4 gas assumptions. For the Mach 10 case, the 
computed C/v values compared well with the tunnel data. However, for the Mach 6 case, there were some 
differences. In the past C at values computed using FELISA had compared well with tunnel data. See for 
example [5]. The differences noticed in the present case could not be explained. 

The general trend of the C 1J} Vs. a curves noticed in the tunnel data are present the computed results. 
The computed pitching moment values for Mach 6 and 10 compared well with the tunnel data at low angles 
of attack. At higher angles of attack the tunnel data show gradual pitch-up compared to the computed 
results at both the Mach numbers. This is likely to be the result of flow separation on the upper surface of 
the canted fin. Flow separation occurs because of the the boundary layer interaction with the trailing edge 
shock. The present computations are inviscid; hence there is no flow separation. However, the trends of the 
tunnel data are predicted by the computations. 
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